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Abstract 

We present a study of both the "Iterative Procedure" and the "/3 — > transfor- 
mation", proposed by Tsallis et al (Physica A261, 534) to find the probabilities pi 
of a system to be in a state with energy 6j, within the framework of a generalized 
statistical mechanics. Using stability and convexity arguments, we argue that the 
iterative procedure does not always provide the right temperature dependence of 
thermodynamic observables. In addition, we show how to get the correct answers 
from the "/3 — > 0' transformation" . Our results provide an evidence that the Tsallis 
statistics with normalized q-expectation values is stable for all ranges of tempera- 
tures. We also show that the cut-off in the computation of probabilities is required 
to achieve the stable solutions. 
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The study of alternatives to the Boltzmann-Gibbs statistics has arisen increas- 
ing interest on the last years. Among several proposals, the most widely used 
and studied has been the so-called "Tsallis statistics" [1]. This new formula- 
tion seems to be appropriate for the study of non-extensive systems [2]. The 
list of non-extensive systems, and that are reported as candidates to obey the 
Tsallis statistics, is growing fast. Therefore, the understanding of the basic 
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principles and properties of this non-extensive statistics has gain fundamental 
importance. In this direction a recent work by Tsallis, Mendes and Plastino [3] 
covers the role of constraints within the Tsallis statistics. In that work, they 
study three different alternatives for the internal energy constraint. The first 
two choices correspond to the ones which have been applied to many different 
systems in the last years [2]. They are: (ia) Y^iPi^i = U and (ib) J2iPi e i — U q . 
However, both constraints present undesirable difficulties. A third choice for 
the internal energy constraint is considered by them [3], 



U q = 



J2j=l Pi e i 

l^i=i Pi 



where q comes from the entropy definition [1] 



Sq — k- 



(2) 



with J2iPi = 1, where i is a given state with energy from fl possible states. 
Each constraint (ia),(ib) and eq.(l) determines a different set of probabilities 
Pi for each state with energy e«. The extremization of the generalized entropy 
(2), under constraint 1 gives us an implicit equation for the probabilities pf 



Pi 
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with 



n 

E 



(i - g)p(6i - u q ) 
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(4) 



The normalized q-expectation value of an observable is therefore defined as 



O q = (Oi) q = 



V s2 r) q 
2^i=i Pi 



(5) 



where O is any observable which commutes with the Hamiltonian - otherwise 
we should use make use of the density operator p. We will refer to this refor- 
mulation of the Tsallis statistics as "with normalized q-expectation values" . 

In order to solve equation 3 Tsallis et al. suggest two different approaches, 
namely the Iterative Procedure and the (5 — > f3' transformation. In the itera- 
tive procedure, we start with an initial set of probabilities and iterate them 
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self-consistently until the desirable precision is reached. In the /3 — > f3' trans- 
formation the set of equations above is transformed to: 



Pl =[l-(l-q)P'e t ]— /Z' 



Z^^[l-(l-g)/3'e,]- 
i=i 



(6) 
(7) 



with 



g 

[l-q)PU q + EUP q j 



(8) 



In order to obtain p i} we go through the following steps: 

(1) Compute the quantities = (1 — (1 — q)j3'ei), Wi G Q. 

(2) If Vi < them Vi = 0; 

(3) Compute ^=Ef =1 l/- /(1 " 9) ; 

(4) Compute Pt ({3') = yl /(1 ~ q) /Z q r, 

(5) Obtain U q (f3') and any other thermodynamical quantities using equation 

(6) Obtain from equation (8). 

This recipe allows the determination of Pi(f3) for all and consequently 

U q ((3) (and any other observable). The second step in the above procedure is 
the well known (and useful) cut-off [3] associated to "vanishing probabilities". 
This cut-off is needed only for q < 1. Because the cut-off is applied before the 
computation of the probabilities, the norm constraint is still respected. 

In ref . [3] , the authors illustrate their ideas by applying them to a system with 
discrete spectrum e n = ne with e > and n = 0, 1,2, ...,N. For N = 1 we 
have the case of two non-degenerate levels. The limit iV — > oo corresponds 
to the quantum harmonic oscillator. Fig. 3 in that paper presents the results 
for the internal energy u q = U q /e, as a function of the temperature, for both 
procedures. The f3 — > f3' transformation generates a non-physical reentrant 
behavior whereas the iterative procedure produces a discontinuity in the free 
energy. These results suggest that the system might be unstable. Therefore, we 
need to figure out the relation between both approaches. We show, through 
a simple example, how to get rid of the reentrant region and consequently 
that the Tsallis statistics with q-normalized expectation values is stable for 
all T. This evidence is specially important because theoretical arguments for 
stability similar to the ones previously used [4,5] are very difficult (if at all 
possible) to be obtained from the implicit equation (3). 
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We consider the same system as Tsallis et ai: a system with discrete spectrum 
e n = ne with e > and n — 0, 1, 2, N. We focus on q < 1 since for q > 1 no 
stability problems occur. We start by computing the generalized free energy, 
F q = U q — ^ nig where ln g (a;) = (a; 1 " 9 — 1)/(1 — g), for a two level system. 
In figure (1) we present the free energy / = F q /e, as obtained from both 
procedures, as a function of the normalized temperature t = T/e. 
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Fig. 1. Free energy calculated for the two-level system with q = 0.3. The curves are 
slightly displaced up only for the purpose of better visualization. The results from 
the (3 — > P' transformation (solid line) show us an unphysical "loop" . There are two 
possible paths for the iterative procedure: the dashed line is obtained by starting 
from a t = configuration. If the initial configuration is the one at t = oo, the 
system follows the doted line. Outside the region shown in this figure, both results 
give the same results. 

Figure (1) shows us that the iterative procedure gives a free energy curve with 
a jump. The jumps will appear at different temperatures if we take different 
initial conditions due to a feature of the iterative method of solving non-linear 
equations: the system tends to go to the nearest solution. The states that are 
reached for only one of the two paths shown in figure (1) should correspond to 
metastable states. We should consider only the one with the lowest free energy. 
In this way, we get rid of the uncomfortable jumps: the free energy cannot be 
discontinuous since its convexity must be guaranteed for all T [6]. Because 
of the simplicity of the system treated here (an energy spectrum with only 
two levels), it is easy to find both paths. However, as we are going to show 
in this work, it is even easier to get the correct behavior from the (5 — > (5' 
transformation. 



The curve for the free energy from the /3 — > f3' transformation displays a 
closed loop. This loop deserves also a careful analysis. The curved line must 
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be discarded by convexity arguments (the specific heat would be negative 
[6]). The remaining states are the metastable states (the same states that 
were obtained from the iterative procedure). Nevertheless, we can also get 
rid of the loop by choosing always the states with the lowest energy. Similar 
procedures and discussions can be found in [7] and [8]. In figure 2, we show 
a closer look on the loops, for some values of q. For this system, we find no 
reentrant behavior for q > 0.56. 
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Fig. 2. Closer look on the loops on the free energy vs. temperature. The curves 
from the right to the left corresponds to q = 0.558,0.5595,0.55975,0.56. The loop 
disappears as q approaches 0.56. 

The procedure we have proposed (to look first at the free energy curve) gives 
the correct behavior for the internal energy and consequently for the com- 
putation of the other thermodynamical quantities besides explaining why the 
iterative procedure and the (5 — > f3' predict different behaviors for the internal 
energy. In fig. 3 we show the corrected curves for the internal energy as func- 
tion of the temperature. For the sake of clarity, we compare this new recipe 
with the j3 — > (3' transformation. For q > 0.56 we do not need to apply any 
further correction because there are no loops on the generalized free energy. 
As shown by Tsallis et al. [3], for those values of q both procedures give the 
same results. The correct expression of the temperature transformation for a 
N = 8-level system is shown in fig. 4. We decided to show this result for a 
richer system then for the simple two-level system, in order to emphasize the 
role of the spurious states introduced by the (5 — > (3' transformation. 

Let us address now the role of the cut-off in the stability of the statistics. 
In fig. 5 we present the free energy calculated with and without the cut-off 
for q = 0.5. These calculations are only possible for q — 1 — l/2n where 
n — 1, 2, 3... since the exponent 1/(1 — q) will be even. As defined in the step 
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Fig. 3. Internal energy as function of the normalized temperature. The thick solid 
line correspond to q = 0.45 internal energy obtained from a (3 — ► f3' transformation. 
This curve should be corrected by the thick dashed line. The dotted line correspond 
to q = 0.5, corrected by the bold dotted line. The dotted-dashed line correspond to 
q = 0.55975, where a reentrant behavior for a small range of t can still be noticed. 
The thin line corresponds to results at q = 0.56, where we do not need to correct 
the /?—►/?' transformation. 
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Fig. 4. t vs t' for N = 8 and q = 0.3. The dashed-line corresponds to the cutting 
points obtained from free energy arguments. We see that t is an increasing function 
of t'. 

(1) of the j3 — > f3' recipe, jji might assume negative values but because they 
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are raised to an even power, y*^ 1- ^ will be always real. Without the cut-off 
step, the free energy shows an unstable concave region. When the cut-off is 
applied, the unstable region disappears. It is clear that the unstable loop does 
not appear due to the cut-off. For these special values of q we were able to 
calculate the roots of eq. (3). The reentrant region corresponds to three real 
and different roots. Outside the reentrant region, we obtain either one real 
and two complex conjugated roots or a three-fold degenerate real ones. 
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Fig. 5. Free energy for the system with (dashed lines) and without (continuous line) 
cut-off (see text). Besides the loop, the solution without cut-off presents a concave 
region wich does not fill the stability requirements. The curves are slightly displaced 
up again for clarity 



In summary, we have shown that the iterative procedure gives a non physi- 
cal discontinuity in the free energy. An alternative procedure called u {3 — > (3' 
transformation" must include a more careful treatment based on free energy 
convexity arguments. We show that the standard way of choosing states with 
lowest free energy values allows us to get a correct internal energy tempera- 
ture dependency. Since the reentrant behavior appears from j3 as function of 
j3', the present prescription restores the proper behavior in the temperature 
dependence of the thermodynamic observable. But the most important con- 
sequence is that the Tsallis statistics with normalized q-expectation values is 
stable, in the sense that, for all /3, the requirements of convexity (concavity) 
for q < 1 (q > 1) in the free energy are satisfied. The current method has been 
also applied to the two-dimensional Ising model in a non-extensive regime[9]. 
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